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THE ADMITTANCE OF A PARALLEL- PLATE APERTURE ILLUMINATING 
A DISPLACED DIELECTRIC-DIELECTRIC BOUNDARY 

By J. Earl Jones 
Langley Research Center 

SUMMARY 

The variational method is used to determine the admittance of a transverse- 
electromagnetic -mode excited, parallel -plate aperture having an infinite flange and illu- 
minating a displaced dielectric -dielectric boundary. The results are applied to provide 
insight into two problems of interest: (a) the effect of an air gap on the admittance of 
dielectric -coated aperture antennas and (b) the admittance of an aperture antenna coated 
with an underdense plasma slab. Numerical computations for the air-gap problem are 
verified experimentally. 

For either problem, the most significant changes in the admittance occur for dis- 
placements less than 0.25 wavelength. (The wavelength referred to is the free-space 
wavelength divided by the square root of the relative permittivity of the medium between 
the boundary and the aperture.) For the dielectric-coated slot antenna, changes in the 
admittance due to an air gap are enhanced if either the slot width is reduced or the per- 
mittivity is increased. For a slot antenna coated with a thin underdense plasma slab, 
the admittance becomes more sensitive to changes in the layer thickness as the plasma 
becomes more dense. As observed theoretically for both the air-gap problem and the 
underdense -plasma problem, the admittance is not influenced by the presence of surface 
waves. Experimental verification of this result is given for the air-gap problem. 

INTRODUCTION 

Flush-mounted aperture antennas covered by low-loss dielectrics are frequently 
used as radomes or for thermal protection of spacecraft antennas on reentry vehicles. 

Air gaps between the aperture surface and the dielectric often result from construction 
errors and unequal thermal expansion of the dielectric and metal surfaces. Such gaps 
cause appreciable changes in the admittance from design values. Other than experimental 
or empirical data, little quantitative information on such effects is presently available. 

One useful model for gaining insight into the effects of an air gap on the admittance 
of a dielectric-coated aperture antenna is the ground- plane -mounted aperture illuminating 



min 


a dielectric layer of finite thickness and displaced a distance d. A two-dimensional ver- 
sion of this model has been treated by Burnside (ref, 1) using the wedge-diffraction method 
from reference 2 for the parallel -plate aperture. Although both finite and infinite layers 
were treated by this method, the results are not applicable for displacements less than 
0.5 wavelength whenever the slot width is less than 1 wavelength (ref. 1) and, hence, are 
not applicable for air gaps of most sizes which occur in practice. 

An alternate approach for obtaining results which are valid for displacements less 
than 0.5 wavelength is the variational method used in conjunction with Fourier transform 
theory (refs. 3 to 8), To achieve some basic insight into the effect of air gaps, the layer 
thickness will be assumed to be infinite in this paper. However, if in the theoretical deri- 
vation both the air layer and the infinite dielectric medium are assumed to be homogeneous 
dielectrics, the results can also be used to determine the admittance of a lossless, under- 
dense plasma-coated aperture radiating into free space; the overdense plasma-coated 
aperture has been studied previously by Fante (ref, 9). 

For both the air-gap problem and the underdense-plasina problem, numerical com- 
putations which show the influence of the permittivity, the slot width, and the layer thick- 
ness on the admittance are given. For the air-gap problem, the results are compared 
with the admittance of the same aperture radiating into a half-space having the same per- 
mittivity as that of the dielectric and also radiating into free space. Similarly, for the 
underdens e- plasma problem, the results are compared with the admittance of the same 
aperture radiating into a half -space having the same permittivity as the plasma and also 
radiating into free space. Experimental verification of the theory for the air-gap problem 
is also given. 


SYMBOLS 


A = 27 ja^ 

a aperture width 

a^ aperture width, free- space wavelengths 

B aperture susceptance 

b normalized aperture susceptance 

Ci^ i ^(k x ) J C2^Hkx) spectral constants 

C{£0 normalized function 


2 


D - 2iTdx 


Da(/3),Dq( 0) normalized functions 

d thickness of medium 1 

d^ thickness of medium 1, free-space wavelengths 

d Xe = /^r,i € r f i d X 
E electric field intensity 

E Fourier transform of electric field intensity 

F(x) function of x 

f{(3) normalized function 

f frequency 

G aperture conductance 

g normalized aperture conductance 

H magnetic field intensity 

H Fourier transform of magnetic field intensity 

H 0 ^(x) Hankei function of the second kind of order zero and argument x 

I reaction integral 

h’h’hM’h integrals 

Im( ) imaginary part of complex variable 

i medium, 1 or 2 {see fig. 1) 

J^^XjQ) surface electric current density on ground plane 
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i = \Fi 


ki 

k x 

k z 

k 0 

Ni 


n 

R 

Re( ) 


SiO) 

t 

u(k x ) 

V 0 


V 

X 

X A 




Y 

Y 


Y 0 


propagation constant in medium i 

x- directed propagation constant 

z -directed propagation constant in medium i 

free-space propagation constant 

refractive index in medium i 

variable of summation 

reflection coefficient (current) 

real part of complex variable 

normalized function for medium 1 

time 

spectral function 
aperture voltage 
complex variable 
dummy variable 
variable of integration 

distance along X-, Y-, and Z-axis, respectively 
aperture admittance 
Fourier transform of Y 
free-space characteristic admittance 
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y 

js 



c i 


normalized aperture admittance 
normalized variable of integration 
propagation constant 

Kronecker delta; 1 when n # 0, 0 when n = 0 
permittivity of medium i 


4 = Re(q) 

4 = Im(ei) 

c r relative permittivity 

e r ^ relative permittivity of medium i 

4,i = Re (er,i) 

e r,i = Im ( e r f i) 

€q free -space permittivity 

= - Nj2 

Aq free-space wavelength 

permeability of medium i 
Mr ? i relative permeability of medium i 

JUq free-space permeability 

4>i = /Ni 2 - 13 2 


a) angular frequency, radians per second 

ojp plasma frequency, radians per second 


5 



Subscripts: 


i medium, 1 or 2 

x, y,z component in and z-direction, respectively 

Superscripts : 

i medium, 1 or 2 

w waveguide 

Abbreviations : 

SWR standing-wave ratio 

TE transverse electric 

TEM transverse electromagnetic 

TM transverse magnetic 


SYNOPSIS OF THEORY 

The geometry of the problem under consideration is shown in figure 1. A synopsis 
of the theory necessary for determining the admittance of the aperture is given in this 
section; mathematical details are given in appendix A. The admittance is determined 
from the stationary expression 

Y = G + jB - r /2 E x ( w) (x,0)H y < w >(x,0)dx (1) 

|v o j -a/2 

where the aperture electric field E x ^ w ^{x,0) just inside the waveguide is given by 
and where the time dependence for the fields is implicit. The boundary -value prob- 

lem for the fields in the external medium is solved with the aid of Fourier transform 
theory* Application of the convolution theorem and the appropriate boundary conditions 
in the aperture plane leads to the result 
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( 2 ) 


where 

and 



ki 2 = 


In appendix B this result is checked for several limiting cases of interest. 

The integral in equation (2) is a contour integral along the entire positive k x -axis in 
the v-plane, where k x = Re(v). Consequently, in evaluating the integral it is necessary to 
include the residue contributions of the real axis poles which physically represent the 
propagation constants of TM (to the X-axis) surface wave modes in the layer. However, 
if the dielectric-dielectric boundary is a "rare-to-dense" type of boundary, that is, if 
eg > e'j for p 2 = 4 1 = 40 and e i = Cjj = 0, no surface wave poles occur and thus no 
energy is coupled into surface waves. Under this condition the integral in equation (2) 
can be easily evaluated by straightforward numerical integration. Also, to facilitate 
numerical computations, the admittance, the constitutive parameters, and the dimensions 
can be normalized to those of free space by letting k x = kg/3. This normalization has 
been done in appendix C, and the final result is 



(3) 


(4) 

b % 2 r ;V i4+: *) 

(5) 


(6) 

h - e r,l*r,2 J Ni D B 0) f( ^ 

(7) 
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(8) 



Dfe 


( e r,l 0 2 J 


cos 


sin 



sin 2D0 1 


f(/3) <*£ 



o sinh Do * 

D (4,i^r° sh do , 

D b {0) 


( e r, 2 ) 2 y- 8inh 2D0 1 


f(0) d(3 


and where 



3j e r,l tanh De i 
+ e^. tanh D0j 


f( 0 ) d 0 


d a( 0) *= (er,1^2) 2cos2D ^l + ( e r,2^l) 2sin2r> ^l 
D b { 0) = (4,1^2) 2cosh2De l + (4,2 e l) 2sinh2Dfl l 


m = 

0 2 

<*>i = \/Ni 2 - 0 2 
0 i = / 0 2 “ Ni 2 
= 

D = 2ird x 

dx “ d/x 0 

a x = a /*o 

0 r ,i - Mi/Mo 


( 9 ) 


( 10 ) 


( 11 ) 

( 12 ) 
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The reflection coefficient (current) is found from the relation 

„ y - 1 
R = 

y + 1 

numerical computations and experiment 

The normalized equations given in the preceding section were programed on a digital 
computer and computations of y and R were obtained for the parameters given in the 
following table: 


Problem , 

Fixed parameters 

Variable parameters j 

Experimental 


Mr,l = M r , 2 = 1 

a x = 0.03 to 0.6 

a x = 0.332 

Air gap 

e r,l ” 1 

c r,2 = 2.05 to 9 

= 2.55, 2.57 


e r,l = c r,2 - 0 

dx = 0 to 2.5 

dx = 0 to 1.0 

Underdens e 
plasma 

Mr,l a M r , 2 = 1 
e r,2 = 1 

ir ft a 

c r,l = € r,2 = 0 

a x = 0.03 to 0.6 ! 

e r>i = 0.05 to 0.8 
dxe = 0 to 1.0 

None 


Experimental verification of the theory for the air-gap problem was obtained with 
the apparatus shown in figure 2. The infinitely thick coating in the theoretical model was 
simulated by a dielectric slab in the experimental model. Two sets of data were obtained, 
one with a plexiglass slab having a thickness of 0.983 in. (2.497 cm) and a relative per- 
mittivity of 2.57 and one with a polystyrene slab having a thickness of 2.010 in. (5.105 cm) 
and a relative permittivity of 2.55. The plexiglass slab was 4.75 wavelengths thick and 
the polystyrene slab was 9.71 wavelengths thick, with respect to the constitutive parame- 
ters of the dielectric at 35.7 GHz, the operating frequency. For either set of data, as 
indicated in figure 3, excellent agreement is obtained in the near field (that is, for 
d^ < 0.5) of the aperture, although the theoretical results include only the dominant mode 
as the trial field. Similar observations regarding the usefulness of the dominant-mode 
assumption in determining the admittance of aperture antennas have also been made for 
both the rectangular (ref. 5) and circular (ref. 6) apertures in the presence of dielectrics. 
The minor disagreement between theory and experiment for dx > 0.5 is believed to be 
caused by the presence of the back wall of the slab in the experimental model (ref. 1). 



Numerical computations of y for the same aperture radiating into a homogeneous 
half- space (ref. 3) for various relative permittivities as a function of a^ (denoted by 
tick marks) are shown in figure 4 and serve as a useful check case for the limiting values 
of d^ e . For dxe = 0 (that is, in the absence of an air gap or an underdense plasma), 
the value of y for the air-gap problem is on the locus corresponding to the relative per- 
mittivity of the coating and at the given value of and the value of y for the 
underdense-plasma problem is on the free-space locus and at the given value of a^. 
Similarly, for d^ e = 00 (that is, an infinitely thick air gap or underdense plasma), the 
value of y for the air-gap problem is on the free-space locus and the value of y for 
the underdense-plasma problem is on the locus corresponding to the relative permittivity 
of the plasma. 

Numerical computations are given in figures 5 and 6 for the range of parameters 
listed in the table. Figure 5 shows the effect of the air-gap thickness d^ upon the 
admittance for various aperture sizes and dielectric constants. The values chosen for 
the dielectric constants in figures 5(a), 5(b), and 5(c), respectively, correspond approxi- 
mately to those of polytetrafluoroethylene (teflon), fused quartz, and alumina. 

In figure 6, the admittance with an underdense plasma slab is plotted as a function 
of the plasma thickness for various combinations of plasma densities and aperture sizes. 
The values chosen for the relative permittivities in figures 6(a), 6(b), and 6(c), respec- 
tively, correspond to plasma electron densities of 0.5, 0.8, and 0.95 times critical density. 

DISCUSSION 

It is evident from the admittance loci in both figures 5 and 6 that the most significant 
changes in y occur whenever the displaced boundary lies in the near field of the aper- 
ture; more specifically, for d^g < 0.25. It is also evident that this result is independent 
of and e^.,2 f° r the air-gap problem (fig. 5) or of a^ and for the underdense- 

plasma problem (fig. 6). For dxe >0,25, the locus, and hence |R(, merely oscillates 
about the value of y obtained from figure 4 for e r = 

For the air-gap problem the amplitude of the |R| oscillations increases with an 
increase in either the slot width or the permittivity, or both. This observation can be 
deduced physically: an increase in ax increases both the directivity and the effective 
aperture so that more energy is collected when reflected from the boundary; also, an 
increase in ej.,2 causes more energy to be reflected at the boundary. The air-gap 
tolerance effect (that is, the sensitivity of y with respect to d^ in the near field) is 
also enhanced for increased values of the relative permittivity but is primarily dependent 
on the slot width, becoming extremely sensitive for narrow slots. As a numerical 
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example, for a narrow slot (a^ = 0.1) at a frequency of 35.7 GHz and a hard dielectric 
boundary ^^2 ~ ®)» iRl is decreased from 0,56 to 0.24 for an air gap of only 0.008 in. 
(0.02032 cm). 

For the unde rdense -plasma problem, as for the air-gap problem, the oscillations 
in |r| as a function of the layer thickness are increased with an increase in the slot 
width for any given value of relative permittivity. Moreover, for a given aperture the 
sensitivity of y becomes more enhanced as the plasma becomes more dense; that is, 
as ojp approaches w or as j approaches zero. 

The admittance locus for = 0.6 and e). j = 0.05 in figure 6 is an interesting 
indication of how a thin plasma slab, close to critical density (^p/^ = 0.975^, will severely 
mismatch a spacecraft antenna previously matched to free space (d^ e = 0). For example, 
}R| is increased from 0.11 to 0.86 for a layer only 1/20 wavelength thick, with respect 
to the plasma parameters. 

An observation of practical Interest concerns the absence of surface wave effects in 
equation (2) for the air-gap problem and the significance of the experimental data in fig- 
ure 3. In a similar problem where medium 2 is a perfect conductor (ref. 8), it was noted 
that the value of y was determined by surface waves. Where medium 2 is a perfect 
conductor, the surface wave modes set up in medium 1 are actually waveguide modes. 

When experimental verification of the theory was sought by use of finite metal plates to 
simulate the infinite plates in the theoretical model, large scatter occurred in the data of 
|R| as a function of d^. The scatter was due to the onset of propagation and the resulting 
mismatch of the modes at the terminations of the metal plates (see ref. 8 for a more com- 
plete discussion) and was evident to some extent even when the edges were curled to 
improve the match. (See figs. 5 and 8 of ref. 8.) It is evident in figure 3, however, that 
no appreciable scatter occurred, although the edges at the terminations were much more 
abrupt than in the metal-sheet problem. These observations further reveal that scatter 
in the experimental data for |R( as a function of d^, when finite surfaces are used to 
simulate infinite surfaces, is caused by the mismatch of surface wave modes at the 
terminations . 


CONCLUDING REMARKS 

A variational method has been used to determine the admittance of a transverse- 
electromagnetic-mode excited, parallel-plate aperture having an infinite flange and illu- 
minating a displaced dielectric-dielectric boundary. The results were used to provide 
insight into two problems of interest: (a) the effect of an air gap on the admittance of 
dielectric -coated aperture antennas and (b) the admittance of an aperture antenna coated 
with an underdense plasma slab. Numerical computations for the air-gap problem were 
verified experimentally. 
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The most significant changes in the admittance were found to occur for displace- 
ments less than 0.25 wavelength. (The wavelength referred to is the free-space wave- 
length divided by the square root of the relative permittivity of the medium between the 
boundary and the aperture.) For the dielectric -coated slot antenna, perturbations in 
the admittance due to an air gap were found to be enhanced if either the slot width is 
reduced or the permittivity is increased. For a slot antenna coated with a thin under - 
dense plasma slab, the admittance was found to become more sensitive with respect to 
the layer thickness as the plasma becomes more dense. 

It was deduced theoretically for both the air-gap problem and the underdens e~ plasma 
problem that the admittance is not influenced by the presence of surface waves. This 
deduction was verified experimentally for the air-gap problem. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., November 19, 1968, 

125-22-02-02-23. 
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APPENDIX A 


DERIVATION OF THE ADMITTANCE EXPRESSION 


The purpose of this appendix is to present the mathematical details leading to the 
development of equation (2). The general geometry of the problem, shown in figure 1, con- 
sists of a TEM mode excited, parallel-plate aperture having an infinite flange and radiating 
into a half -space environment consisting of (1) a lossy homogeneous dielectric layer of 
thickness d and with constitutive parameters p j and ej = e'j - jej and (2) a lossy 
homogeneous infinite medium with constitutive parameters p 2 ^d «2 = e 2 “ i € 2* 
time convention for the fields is used throughout the deviation. 

An expression for the admittance of a dominant-mode excited, ground-plane- 
mounted, cylindrical waveguide of arbitrary cross section is derived in reference 3 
(pp. 8-12) in terms of the aperture electric and magnetic fields just inside the waveguide. 

It is further shown (ref. 3, pp. 12-16) that this expression is stationary with respect to 
small variations in the true value of the aperture electric field, provided the fields in the 
environment into which the aperture radiates satisfy the radiation conditions. For the 
aperture in figure 1, this expression reduces to 


Y 



(Al) 


where I, the reaction integral, is 


I 



E x ( w >(x,0)H y (w) (x,0)dx 


(A2) 


and where the aperture field is assumed to consist of the incident and 
mode; that is 

E x < w >(x,0) (lx I 

E x (w) (x,0) = 0 (|x | 

The normalizing factor -2— in equation (A3) is necessary in order that the orthonormal 

vEL 

property for the TEM mode function be satisfied. (See, for example, ref. 10.) 

The reaction integral in equation (A2) can be evaluated by considering the assumed 
aperture electric field to be an equivalent magnetic surface current source for the fields 
in the external environment (z > 0) and solving for the magnetic field arising in the aper- 
ture as a result. For the two-dimensional problem under consideration, it is well known 


the reflected TEM 


s i) 


(A3) 
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that the Maxwell equations can be decoupled into two independent sets, one TE and one TM 
to the Y axis. Thus, with TEM mode excitation, the fields in both mediums 1 and 2 must 
satisfy the TM~decoupled wave equation 

f>V‘W> + ^h/W) 2 (i) 

a z 2 1 * 

and the TM-decoupled Maxwell equation 


9x 2 


+ kj 2 Hy^'jXjZ) = 0 


ft) 


where 


E x v ; {x,z) = - ri— 
3&*i 


j 8H y ( i )(x,z) 


Bz 


and i = 1 or 2. 

Equations (A4) and (A5) must also satisfy the following boundary conditions: 


(A4) 

(A5) 

(A6) 


E x (1) (x,d) = E x < 2 >(x,d) 


(A7) 

H y <%,d) = H y < 2 >(x,d) 


(A8) 

E X (1 W) = E x (w) (x,0) = ^2 

,/a 

"SicsT 

VII 

(A9a) 

E x ^>(x,0) = 0 

( lxl> f) 

(A9b) 

H y <%,0) = Hy^ w )(x,0) 

( lxl 5 1) 

(AlOa) 

H y ^(x,0) = -J x ( x )(x,0) 

( ,xi> !) 

(AlOb) 


as well as the boundary condition requiring only outgoing waves at infinity. All these 
boundary conditions are readily enforceable except equation (AlOb), since the ground- 
plane surface current arising as a result of the source field in the aperture cannot readily 
be determined. An alternate boundary condition is achieved by requiring the conservation 
of reaction (ref. 11, pp. 420-421} across the aperture. Thus 


I = r a//2 e x ^ W \x,Q)H ( w ^(x,0)dx = P a / 2 E x ^)(x,0)H v ^{x,0)dx (All) 

J -a/2 J -a/2 J 


The boundary -value problem for the fields in the external medium is easily solved 
by assuming Fourier transform solutions of the form 
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H y (i >(x,z) = ^ §°° H y (l) {k x ,z)e“ ;ik x x dk x 


(A12) 


thus reducing the number of independent variables in the wave equation from two to one, 
yielding the total differential equation 


where 


dW) . 2iI (i), 

‘ ^2 H y ( k x>zJ - 


k ,2 _ t.2 _ u„2 

K z,i - K i K x 


0 


(A 13) 
(A 14) 


Similarly, substitution of equation (A12) into equation (A 5) gives the transformed field 
equation 


E x ^( k x,z) = - 


1 dHy (i) (k x> z) 
j dz 


(A15) 


and substitution of equation (A 12) into equations (A 7), (A8), and (A9) yields the transformed 
boundary conditions 


Ex (1) (k x ,d) = E x { 2 } (k x ,d) 

(A 16) 

H y < 1 \k x ,d) = H y ( 2 >(k x ,d) 

(A 17) 

Ex (l) (k x ,0) = E x (w)( kxj0 ) 

(A18) 

The general solution of equation (A13) is 


H y (i >(k x ,z) = Ci(i)(k x )e Jk =’ iZ + C 2 ( i >(k x )e^ kz * iZ 

(A 19) 

which yields, when substituted into equation (A 15), the result 


E x (i) (k x ,z) = - ^ Ci (i) (k x )e ]kz A Z + ^ C 2 % x )e‘ jkz > iZ 

(A20) 

where the branch 

k z,i = ( k i 2 - k x 2 ) (Re(k Z) i) > 0; Im(k Zji ) £ 

O) (A21) 

must be chosen to satisfy the radiation condition for large values of z. 

The spectral constants Ci^)(k x ), C 2 (l)(k x ), Ci(2)^ kx ), and C 2 (2)(k x ) are eval- 
uated by applying equations (A16), A(17), and (A18) and requiring only outgoing waves for 


large z. This latter condition requires that Cj(^^k x ) = 0. Using this result and 
applying equations (A16) and (A17) to equations (A20) and (A19), respectively, give 
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where 


- Ci a>(k x )e )k M d + c 2 (^k x )e*.> d . € -i^ic 2 ( 2 )(k x) e- )k2 . 2d 

€2 k 2,l 

Ci^^(kx)e^ z >^ + C2^\k x )e = C2^ 2 ^kx) e 

Eliminating C 2 ^(k x ) from equations (A22) and (A23) gives 

C 2 (1) (k x ) = U(k x )Ci< 1 >(k x ) 

, ^2,2^1 

, , 1 + j 2 k z id 

U(k x ) = tli_2 e ’ 

**7 


1 - 


C Z, 26 l 

kz,l e 2 


(A22) 


(A23) 


(A24) 


(A25) 


Substitution of equation (A24) into equations (A19) and (A20) yields 


and 


From the definition 


it is evident that 


= [>M 2 + U(k x )e-J k M^C 1 <D(k x ) 

(A26) 


(A27) 

HyW^z) 
*(“*>*) ' _ v 

(A28) 

u*e, > 12 - U(k x) a- Jk ^ 12 

1 ^ 

(A29) 


The solution is completed by applying the boundary conditions given in equations (A9) 
and (All) in the Fourier transform domain. From equation (A9b) the limits of integration 
on the right side of equation (All) can be extended to infinity without affecting the result. 
Furthermore, because of symmetrical excitation of the aperture, the fields everywhere 
are even functions of x. Then, by the convolution theorem 
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Substituting equation (A28) into equation (A30) gives 



1 = h I* Y ( kx >°) (% <1) ( k x»°]l 2 dk x 

(A31) 

But by virtue 

of equation (A9) 



E x (1) (k x ,0) = E x < w >(k x ,0) 

(A3 2) 

and 

k x a 




(A33) 


Consequently, from equations (Al), (A31), (A32), and (A33), the admittance is found to be 


y -M" y <**- o > 


sm 2 

k * 2 


dk* 


(A34) 


Equation (A34) is a particularly important result , since the admittance of the parallel - 
plate guide is now expressed as an integral over mode space and involving the spectral 
constant Y(k x ,0), which is dependent only on the nature of the external environment. 
Furthermore, since both Y(k x ,o) and its inverse transform 


Y(x,0) = 1- Y(k x ,0)e" ]kxX dk Ji 


(A3 5) 


are even functions of k x and x, respectively, the integral in equation (A34) can be con- 
verted from mode space to physical space by use of the convolution theorem to give 


Y-if Y(x,0) F(x) dx 


where 




.00 sin 


k*a 


2 22 


k x 2 


2 “ J k x x 

” e dk x 


(A3 6) 


(A37) 


Application of residue theory to equation (A3 7) yields 

F(x) = I (a ~ |x|) 

F(x) = 0 

Thus, an alternate form for equation (A34) is 

Y = | P (a - x)Y(x)dx 

a J 0 



(A38) 


(A39) 
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Equation (A39) may be more useful for purposes of numerical computations since the range 
of integration is finite. 

For the problem under consideration 


or 


Y{k x ,0) = 


cuej U{k x ) + 1 
1 


Y(k x ,0) = 




A ' kz > 2e l 
f 1 +] k z,P 2 


C M 


‘ k z,l € 2 



(A40) 


(A41) 


The inverse transform of equation (A41) is not readily evaluated except for certain special 
cases; consequently, in the spatial frequency domain the admittance is given by 


Y 




sin 2 5^ 

L + j ^ 7l ,2( 1 tan k z id\ 
wei / k z> ie 2 

kx 2 

M^ + jtan * J 
V%P2 / 


dk 


■x 


(A42) 


which is the result given as equation (2). The surface wave propagation constants are then 
found by solving the equation 

■ Z ’ 2 *■ + j tan k z id = 0 (A43) 

k z,l«2 

for the real roots of k x . 
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APPENDIX B 


LIMITING CASES 


As a partial check on the analysis, equations (A41) and (A42) are compared with the 
following special cases, for which Y has previously been determined: 


(1) Medium 1 a lossy dielectric and medium 2 free space - Under this condition 
equation (A41) becomes 


Y(k x ,Q) 


(t)6i 


1 * J i T7T tan kz ’ ld 
K z,i e o ' 


k Z) 2 e i 


k 


M £ 0 


+ j tan k Zj id 


(Bl) 


where k Zj 2 = (ko 2 - kx 2 ) 1 ^ 2 with kg = W\JWq£o, R e(k Zj 2) > and I m (kz,2) < 0. This 
result is identical with the result established by Compton (ref. 3, pp. 39-40). 

(2) Width d of medium 1 allowed to approach zero - For this case the external 
medium is a lossy half-space with the constitutive parameters pg and £2* Under this 

u>eo 

condition equation (A41) approaches the value and thus equation (A42) becomes 

k Z} 2 


Y 



^ »a 2 
k x 2 k z,2 


dk x 


(B2) 


Equation (B2) is the result established by Compton (ref. 3, p. 25). For this case Y(x,0) 
is found to be (ref. 12, pp, 30-42) 



we 2 -jk x x 

i e 

k z,2 


^x = ^H 0 (2) (k 2 (x|) 


and thus, from equation (A40), the alternate form is 


(B3) 


Y(x,0) = ( a ~ x)H 0 ^ 2 \k2x)dx (B4) 

(3) Dielectric constants of medium I and medium 2 equal — For this case the exter- 
nal environment is a semi-infinite dielectric with the constitutive parameters Pi = P2 
and ej = C2- Furthermore, as ej approaches 63 and M - 1 approaches 11 2 > the bound- 
ary effectively disappears; consequently, the admittance should reduce to the value given 
by equations (B2), (B3), and (B4), independently of d. This condition is evident from 
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examination of equation (A41). Note that, since k 2 \ = k z 2 for this case, the quantity in 
parentheses is equal to unity for any d, and hence 


Y < k *’ 0) “CTid 

as required. 

(4) Medium 2 a perfect conductor - Under this condition the aperture is illuminating 
a perfectly conducting wall. For this case the imaginary part of e 2 is infinite but 

k z ,2 _ (k 2 2 - kx 2 ) 1 / 2 _ (oj 2 p g e 2 - k x 2 ) 1/2 


e 2 


6 2 


e 2 


(B5) 


or 


kz,2 _ |(^ 2 P 2 e 2 - k x 2 ) - ju)2p 2 e|] 

= ^ - J% 


1/2 


<B6) 


tt kz 2 

and as e 2 approaches infinity, the ratio approaches zero. As a result, equa- 


tion (A41) becomes 


and equation (A42) becomes 


e 2 


Yk x ,0 = - 


^ e l 


k Zj i tan k z>1 d 


(B7) 


= . 4 ir 

7raJ 0 


■<* sin 


;^2 k X a 




kx^ t&n k z,l d 


dk^ 


(B8) 


which is in agreement with the result established by Jones and Swift (ref, 7, pp* 13-14), 
who have also converted this result to the alternate representation. For this case (ref. 7, 
pp. 31-35) 


where 


Consequently 


Y(x Ol.lf — 1^1 , - 
v,/ 2jr J_oo k Zjl tank Zjl d 

oil/2 


-jk x x 


™ "jy n fxi 


dk x = we a 


n=0 i 1 + 6 0 n )>n d 


(B9) 




2l v 


? r a v e" )r ” X 

Y=lj o ( s -^ Irv 


dx 


(BIO) 


(Bll) 
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Equation (Bll) is readily simplified to give 

ntoiC- 

Physically, each of the n terms in equation (B12) represents the admittance presented 
to the aperture by the nth TM (to the X-axis) waveguide mode in medium 1; the y n in 
equation (BIO) is the propagation constant for the nth mode. 


“SiS — 3H 1 - cos *n a ) - K y n a - Sin >n a 3 
6 oV 


(B12) 
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APPENDIX C 


NORMALIZATION OF THE ADMITTANCE EXPRESSION 


For convenience in numerical computations, the results given in equations (A42) and 
(A43) can be normalized by letting 

k x = koi9 (Cl) 

Then equation (A42) can be manipulated into the form 


where 


Y = v~ = 
Y 0 


Y G + jB 


Y 0 


= g + jb = 


2e 


r.I p 
'-’O 


77 2 ay 


c(g) 

s x O) 


f (j8) d/3 


(C2) 


. o j3A 
sum 



*<0) = (0 * 0) 



y3 a 

> 

(C3) 

m = {*a x ) 2 If) = 0) 



A = k ft a = —• a = 2irax 
x 0 


<C4) 

, „ , . sin DSi(fl) 

m cos + ,er > 1 “ 2( ' 3) d S iW - 

s l(0) e r iS 2 (0)cos DS x (0) + je rj 2 s i(0) sin Ds/jS) 


(C5) 

D = k Q d = ^ d = 2 Tidy 


(C6) 

1/2 / 

Si</3) = (Ni 2 - 0 2 ) (Re^p] > 0; Im[Si(p 

S °) 

(C7) 

Nj 2 = P r>i 6 r) i 


(C8) 

T tt 

I , Tt £j . e i 

E r,l ' e r,i - )tr,l = ij = ^ ' 3 £ ~ 


(C9) 





(CIO) 


and, in normalized form, the surface wave poles correspond to the real values of /3 satis- 
fying the equation 
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e r oSitjS) 

tanDS 1 0)=j r *Y,l (Cll) 

er,l s 2(/3> 

The results given in equations (C2) to (Cll) can be further simplified and applied to 
the problem of interest in this paper by assuming that both medium 1 and medium 2 are 
lossless = e^2 = 0), nonmagnetic = Pr,2 “ *) dielectrics where = e T j S 1; 

that is, Ng = Nj. An examination of equation (Cll) under these conditions reveals that 
no surface waves are excited in the dielectric layer. Then by use of equation (07) in 
equations (C2) and (C5) the admittance can be written in the form given by equations (3) 
to (12). 

For the aperture radiating into a lossless dielectric half-space, equation (B4) is 
given in normalized form as 


y = 2 ® r C(’ " (C12) 

This result is used to obtain the computations in figure 4. 
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REFLECTION COEFFICIENT MAGNITUDE 


— THEORY (variational solution) 

A EXP. (0.983- in. (2.497-cm) plexiglass slab; £ r = 2.57) 
O EXP. (2.010-in. (5.105-cm) polystyrene slab; e r = 2.55) 

FREE- SPACE LEVEL (theory) 

FREE-SPACE LEVEL (experiment) 



Figure 3.“ Comparison of theory and experiment for reflection coefficient magnitude as a function of width of medium l 


















t fi 2 ~ 2,05. Concluded 
figure 5,- Continued 





(b> e[* 2 ” 4. 
Figure 5.- Continued. 
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(b) £ l Tf 2 - 4* Continued, 
Figure 5 >- Continued. 








(b) tf'2 ~ ^ Concluded, 
Figure 5.- Continued. 
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